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Abstract 

We consider optimally spin-squeezed states that maximize the sensitivity of the Ramsey spec- 
troscopy, and for which the signal to noise ratio scales as the number of particles N. Using the vari- 
ational principle we prove that these states are eigensolutions of the Hamiltonian H (A) = XS^ — S x , 
and that, for large N, the states become equivalent to the quadrature squeezed states of the har- 
monic oscillator. We present numerical results that illustrate the validity of the equivalence. 
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I. INTRODUCTION 



Coherent control of ensembles of two-level atoms may lead to the production of entangled 
spin-squeezed states Q. Each atom of the ensemble can be treated as a spin 1/2, and, in 
analogy with squeezed states of photons P and phonons for squeezed spin states the 
quantum noise in the different components of the total spin is redistributed preserving the 
uncertainty relation. Depending on the context, different definitions of spin squeezing can be 
used. Starting from the uncertainty relation AS x AS y > \(S z )/2\ (and cyclic permutations), 
a possibility is to define states satisfyingA 2 ^ < \(Sj)/2\ as spin-squeezed j^]. Another defi- 
nition was proposed by Wineland et al. |5[ , who showed that the resolution in spectroscopic 
experiments on iV two-level atoms is determined by the factor 

e = ^ (i) 

K \(S)V w 

which measures the quantum noise in a direction perpendicular to the mean value of the 
total spin. States with small £ are relevant to quantum information and constitute the focus 
of our paper. We will also discuss Kitagawa and Ueda's[l| squeezing parameter given by 

£ ■ - (2) 

f.p» l(s)r (2) 

as well as the parameter £ 2 proposed by Raghavan et a/.pj: 

£ = ASl (3) 
?2 S 2 -(S Z ) 2 ' {) 

From the commutation relation for the spin components in three orthogonal directions 
we have that 

(S> AS » * 2" « 



Since the fluctuations for N spins 1/2 are bounded [AS < y/N/2{N/2 + I)], there is a 
lower bound (the Heisenberg limit) for the squeezing factor: 

^4 < 5 > 



for large N. On the other hand, for product states (or unentangled states) £ = 1/y/N. The 
simplest example is the completely polarized state in the x direction | — >, — >, — >, • ■ • , — >), 
which has (S x ) = N/2, (S y ) = (S g ) = 0, and (AS Z ) = (AS y ) = s/N/2. 



Given the above bounds, the key to approaching the Heisenberg limit is to find states 
with | (S x ) | of order N and a fluctuation AS Z of order one. 

Consider for example N = 3 with S — 3/2. The Hilbert space has dimension 4, and we 
can label the states by their projection m = —3/2, —1/2, 1/2, 3/2 along the z direction. The 

Q 

simplest squeezed state 1^/5=3/2) has the form [2J 

|l/2) + |-l/2) 
1*5=3/2) = -j= • (6) 

This state is polarized in the x direction, with (^3/2! £11*3/2) = 1 , {^E f 3/2|»S'y 1^3/2) = 
(*3/ 2 |S , ,|* 3 /2) = and (^3/2|#^3/2> = 1/4, giving £ = 1/2 < 1/^3. 

If we include the states with maximum projection (S z = ±3/2) we will increase both 
AS Z and (S x ) and the question is whether the inclusion of these states can improve the ratio 
AS Z /{S X ). The natural step is to consider a state of the form: 



1*5=3/2) 



|l/2) + |-l/2) |3/2) + |-3/2) 
+ a- 



y/2 V2 



(7) 



for which 



„ 1 J(l + a 2 )\l + (3a) 2 } 1 , / — . , 

This result indicates that the inclusion of larger spin projections with a small amplitude 
can improve the spin squeezing. For a = l/v3 the state has well defined S x = 3/2, and 
£ = l/y/3. For a ~ 0.16 the squeezing factor has its minimum value £ ~ 0.44. 

In section |H] we prove, using the variational principle, that the states that minimize £ 
(the optimally squeezed states) are solutions of the Hamiltonian 

H(X) = XS 2 Z - S x , (9) 

where, without loss of generality, we chose the x axis in the direction of the mean spin and 
S± = S z . 

For different values of the parameter A, the eigenstates of Q have the minimum AS Z 
for a given (S x ) and, in the large iV limit, the squeezing factor scales as £ ~ l/N. For 
1 <^ AS Z <^ VN, H(X) can be mapped to the harmonic oscillator Hamiltonian, and the 
optimally squeezed states become gaussians of reduced variance, identical to the quadrature 
squeezed states. In section IIHI we present numerical diagonalizations of for up to 600 
spins to illustrate the validity of the mapping. 



II. VARIATIONAL CALCULATION 



For N spins, the general state with (S) in the x direction can be written as the following 
superposition of Dicke states \S,m) (eigenstates of the total spin S and the spin projection 
m in the z direction): 



m = J2 ( fs, m \S,m), (10) 

m,S 

with ips,m real to guarantee that (^\S y \^/) = 0. The squeezing factor can be written as the 
following functional of the set {<fs,m}' 



£{Ws,m}) 



1/2 



1/2 



Es m ( PS,m^S,m ( PS,m+X 



(11) 



with m integer (half integer) for N even (odd), and 

X s , m = y/(S + m + l)(S-m) 
twice the matrix element of the operator S x . The minimization condition 

0<fS,m 

gives rise to the following difference equation for (ps,m 



(12) 



(13) 



m—l i 



E(\)<p 



S.mi 



(14) 



with A and E(X) functions of the sums in Eq. ()ll)l . Eq. ()14jl is identical to the eigenvalue 
equation for the Hamiltonian H(\) = XS 2 — S x , and in order to find the minimum squeezing 
we can regard A as a parameter and minimize £ with respect to A for each eigenstate of H(X). 
Also, since H (A) is even in the projection in the z direction, the solutions are even and odd in 
m, and (^IS^Nl/) = 0. It is interesting to note that Eq.(|14|) corresponds to a one dimensional 
tight-binding Hamiltonian for a "particle" in a parabolic potential of strength A and hopping 
to near neighbors given by Xg y m- The squeezing factor can then be interpreted as the ratio 
of a mean square displacement and a mean kinetic term. For large S the parabolic potential 
confines the particle to small values of m and as a result (m 2 ) ~ 0(1), and (S x ) ~ O(S), 
which implies that £ is minimized by states corresponding by the maximum total spin. We 
work in this subspace and in what follows we omit the subindex S in ips, m - 



Hamiltonian (j^J) was considered previously by Law et al. to generate spin-squeezed states 
and also arised several times in the context of Bose-Einstein condensation. Specifically, 
if we write the spin operators in term of the Schwinger bosons a* and b^\ 

S z =- {a) a - b ] b) , S + = alb, (15) 

Hamiltonian Q becomes (ommiting terms that depend on the total number N) 

H(X) = 2A [(a f a) 2 + (fr+fr) 2 ] - l - (a% + b ] a) . (16) 

The Hamiltonian above corresponds two condensates, each of them treated in the single 
mode approximation, with the S x term in @ becoming a Josephson coupling betwen the 
condensates. This correspondence was used by Milburn et a/.j^, Steel and Collett [ill , Gor- 
don and Savage and Raghavan et a/.p]. In this representation, the squeezing parameter 
£ can be interpreted as the ratio of the particle number fluctuation (AS Z ) to the Josephson 
energy {{S x )) of two weakly coupled condensates: start with two condensates (a and b) with 
equal particle number N a = = N/2. This corresponds to the Dicke state \N/2,0) with 
S z = 0. Now we turn on a small Josephson coupling between the condensates that will induce 
a small fluctuation in the particle number of each condensate A(N a — Nb) = AS Z ~ 0(1). 
On the other hand, for the Josephson energy we have 

(ab ] ) ~ (a) (&*) ~ O(N), (17) 

giving £ ~ 1/iV. 

Before analyzing the detailed form of the solutions of ((Hj) we can gain insight on the 
optimally squeezed states by examining t he g eneral behavior of E(\) for small and large A, 



and using the Hellman-Feynman theorem 
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(18) 



dE _ I dH 
~d\ ~ \~dX 

to compute the averages: (m 2 ) = dE(X)/d\ = (A^) 2 , and J2 m ip m X s ^ m ip m+1 = 
XdE(X)/d\ - E(X) = (S x ), giving 

~ dE(\) ~\ 1/2 

£(A) = 1^ J . (19) 

X d -T-E(X) 

For small (m 2 ) (large A), the solutions for even and odd iV are qualitatively different. For 
A — > oo and even, the ground state energy the solution corresponding to the "particle" 



at the bottom of the potential well (m = 0). Using perturbation theory keeping the lowest 
lying states (m = 0, ±1), 



for iV even. On the other hand, for odd N the particle can be either at m = +1/2 or m = 
— 1/2, the energy diverges as A/4, and the ground state corresponds to a linear combination 
as in Eq.(jnj). The condition d^/dX = implies 

§T^=°- ( 21) 

Since E(X) has monotonic curvature, the minimum is given by the product XE(X). Given 
the discussion above, for N even the minimum value of £ corresponds to A — > oo: 

£ (min) = ^, (even N). (22) 

This absolute minimum is however not very interesting since it corresponds to a ratio of 
zeroes therefore implying a vanishing mean value of the spin. This limit has been pointed 
out by Andre and Lukin for a perturbative calculation. For odd N, since E(X) diverges 
for large A, the product XE(X) changes sign and the absolute minimum of £ corresponds to 
a large but finite value of A that we determine numerically by solving exactly Hamiltonian 
©• 

If we extend our variational treatment to the squeezing parameter given by Eq.(J2J), we 
see that the states that minimize £ sp i n are also eigenstates of ©, but this time we have to 
minimize £ sp m(A) given by 

~ dE(\) ~\ 1/2 

£s pin (A) = r } dX j — 2 . (23) 



dX 

The minimum value of C, sp - m also corresponds to A — > oo, but the square root in the 
denominator of (J23|) implies that 

& = °" (24) 
On the other hand, the variational calculation shows that the states that minimize £ z 
given by Eq.® are eigenstates of 

H S (X) = XS 2 Z - S z . (25) 



As discussed in Ref. p, the states that minimize £ 2 are eigenstates of S z . 



III. NUMERICAL RESULTS 

The numerical minimization for different A (see Figure shows that, in the large A 
limit: 

Un = ~jf, (odd A). (26) 

In general, for larger values of £ the solutions of ((H]) give the optimally squeezed states 
with minimum AS Z for a given value of (S x ), and for which £ scales as 1/A. In order to 
show this, we first note that the exact normalized eigenstate of S x with eigenvalue A/2 and 
A even |l5| has the form of a gaussian of width y/~N : 

„^>=2-n "j^JA'V/, (27) 



iV 

2 

f-iV/2) _ / ,N m , W/2) 



The eigenstate corresponding to = —A/2 is simply <y2m = (— l)" 1 ^™ • Using 
the explicit expressions for X^/2,m it is straightforward to show that the states with smaller 
projections S x = N/2 — n are of form ifm^ 2 oc hffl (m)(pm , and hg(m), with 5 = A/2, 

2 

obeying the following recursion relation: 



N 

T 



4 n) (™+ 1) + hP(m- 1) -— /i^(m+l)-/i^(m-l 



m 



>V™ _L -n _ *»/ 



— — n ) h^h \ in i 



(28) 

For large A we can take the continuum limit of (m) and ()28|) becomes the differential 
equation satisfied by the Hermite polynomials H n (x) [16], with x = m y/2/N. This proves 
that, in the large A limit, the eigenstates (fm^ 2 of S x in the S z representation are harmonic 
oscillator wave functions ip n (x), with the following operator equivalences: 

A 



Sr. — > W —X, 

A + 1 1 / rf 2 A 

+ 2(-^ + I J- (29) 

and an implicit truncation in the allowed excitations n of the harmonic oscillator Hilbert 
space: n < A/2. [For quantum numbers A/2 < n < A, the solutions are of the form 



^(N/2-n) _ ^_x) m (/ 9 7V _ n ( miv /2/A^).] The operator equivalence of Eq. (}2T?|) corresponds to the 
well known Holstein-Primakov transformation 17| provided AS X /N -C 1. Also, the corre- 
spondence between the Harmonic oscillator algebra and that of the angular momentum was 
discussed using the so-called group contraction by Arecchi et al. In this formulation, 

(JHJ) becomes the harmonic oscillator Hamiltonian and the ground state is a rescaled gaussian 
<Po(x) = T[~ x l^o~ x l 2 exp —x 2 /2a 2 , which, in the discrete representation has the form 



1 



^ m ~ —e~^ 2AS ^, (30) 

(2vr)V4A^/ 2 1 J 

with AS Z = u^/N /2. Since the mapping to the continuum is exact provided 1 <C AS Z <C N, 
taking the limit of large N with AS Z fixed we obtain 



I/O m N+1 1 ( 2 1 \ N 



8(AS Z ) 2 \ ■ (31) 
In Figure El we present comparisons of the continuum approximation with exact diago- 
nalizations of dHj) that show excellent agreement for large N and AS Z as small as 0.5. In 
Figure El we show comparisons of the exact eigenstates <p m of Eq. (JHJ) with the analytical 
solution. 

Notice that the operator S y in the discrete (or tight-binding) representation has the form 
of a velocity operator: S y = (i/2) J2 m (-^-N/2,m I m + l)(m\ — h.c). In the tight-binding limit 
the commutator of the position (corresponding to S z ) and the velocity operator is not a 
constant but is proportional to the kinetic energy operator (corresponding to S x ) and the 
commutation relations are equivalent. In the continuum approximation S y — > ^y^S and 
AS y is given by 



^ = <--— ) =7T7T, (32) 



~2dx^/ ~AAS Z ' 

which proves that the optimally squeezed states in the large N limit (and for 1 -C AS Z <C N) 
correspond to minimum uncertainty wave-packets satisfying AS z AS y = (S x )/2. Note that 
the optimally squeezed states given by Eq. (|30|) can be written in the form 

%(N,M,a) = C exp (—aS 2 ) exp {—mS y /2)\N/2, M), (33) 

with a = 1/(2AS Z ) 2 - 1/JV ~ 1/(2AS Z ) 2 and M = N/2. The state ® has some similarity 
with the minimum uncertainty states considered by Rashid [3, Q , but the S 2 (as opposed 

8 



to S z for the Rashid states) in the first exponential guarantees that the squeezing is in the 
direction perpendicular to the mean spin for all M. 

We close by noting that the harmonic oscillator equivalence extends to the dynamic 
generation of squeezed states by the application of an impulsive non-linear Hamiltonian. 
Consider an ensemble of identical two-level atoms with energy splitting Huq. We define 
the corresponding spin quantization axis in the x direction so that H Q = u S x , and the 
equations of motion for S z (t) and S y (t) are exactly the same as those of x and p for a harmonic 
oscillator of frequency u> . If a pulse of the form H' = 5(t)r]S^ is applied, the quasiprobability 
distribution in the (S y , S x ) plane is modified as in Figure EJ the only difference with the 
harmonic oscilator being that 77 < l/\/N (for larger rj the response is periodic in rj). For 
times t > the distorted distribution rotates at frequency u and the squeezing factor 
evolves as 

1 [l + (77iVsina;ot) 2 -77iVsin2a;ot] 1/2 

If we call rj = ao/^/N, with a < 1, the minimum squeezing £ min = ([1 — al/2]a N)~ 1 
scales as 1/N and is reached twice during the cycle of the rotation of the ellipse of Figure 0] 
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Appendix 

In this section we will derive the recurrence equation for the discrete polynomials of Eq. 
(J2EJ)- The normalized eigenstate \n) x of S x corresponding to the eigenvalue S — n, with 
S = N/2, can be written as 



n) 



l)x = Vf ) ( m )\ m )zi ( 35 ) 

m=—S 

with \m) z the eigenstate of S z with eigenvalue m. For our choice of quantization axes we 
have 



S x \n) a 



(S+ + S-) 



n) x 

n — 1 
" 2S 



1/2 



n \ 1/2 



u\n -l) x + (1- —) + l\n + 1) 3 



(36) 
(37) 



which, substituted in (J35j) gives 



2( J^l^mj 



1 - 



n — 1 
~2S~ 



1/2 



2n^" 1) (m) + (l - ^) V2 v^TTj^M. 



Now we define hg(m) as 



<Ps (™) 



V S — m ! ' 



V2"r?T 

and subtitute in ()38|) to obtain the following recursion relation 



(38) 



(39) 



4" +1) (m) = 2 / "' 



1 



ha(m) — 2n\ 



'\-{n-l)/2S, n „ x) 



h [ s L> (m) (40) 



.^fs) y/l-n/2S' 3 v " v """V l-n/25 

Note that, in the limit 5 — > oo, m/\/~S — > x, hg (m) — > H n (x), and (jjQj) becomes the 
recurrence equation for the Hermite polynomials: 



H n+1 (x) = 2xHJx) - 2nH n -x(x). 



(41) 



Iterating (jjOj) we obtain the first few discrete polynomials hg(m 
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hg\m 
(m 



hg(m 



m 



Vs 

7s) 



4(^) -2 



A/1-1/25 



V(l-l/25)(l-l/5) 

tI) 4 - 48 



I (j ( ] - 48 f^=V (1 - 2/35) + 12(1 - 1/5) 



(42) 



v/(l-l/25)(l-l/5)(l-3/25) 
Finally we note that the above discrete Hermite polynomials are proportional to the 
Karwtchouk polynomials K n (x;p, N) 21], which obey the following recursion relation: 



p{n-N)K n+1 (x;p,N)-[n(2p- 1) + x - Np]K n {x;p,N)+n(p-l)K n „ 1 {x;p,N) = 0. (43) 



Comparing ()43|) with (j4"Uj) we see that 



h^\m) 



2n/2 



y/(2S-n)\ 



K n {m + 5; 1/2, 25). 



(44) 
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0.852 




FIG. 1: Numerical results of the absolute rninimum £min for even iV (ctrid S — cis &- function 

of N. 





FIG. 2: AS" Z vs (S^) for the optimally squeezed states for even and odd number of spins N obtained 
by numerical diagonalization of Ea. Q14j) . The end points for TV" even and odd correspond to the 
eigenstates of S x (with eigenvalue N/2) and AS Z = \HV /2. For odd N there is a minimum value 
of (S x ) = (N + l)/4 that corresponds to AS Z = 0.5. This is the bifurcation identified by S0rensen 
\l\ for N = 3 and discussed in the text. The dashed lines correspond to the analytical result 

AS Z = (l/y/8)/y/l-2{S x )/N. 
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FIG. 3: Comparison between the exact eigenstates of H = XS Z — S x for 600 spins (circles) and the 
solutions of the discretized harmonic oscillator of Eq. lJMOjl (x's) in the continuum approximation. 
The curves correspond to AS z /AS zfl = 0.03 (a), 0.08 (b), 0.24 (c) and 1.0 (d), with AS Z>0 = VN/2. 
For AS Z > 0.3AS Z fi the two solutions are indistinguishable. 




FIG. 4: Quasiprobability distribution in the (S z ,S y ) plane for iV spins, before and after a pulse 
H' = 5(t)r]S z is applied on the lowest eigenstate of Hq = ujqS x . In contrast with the harmonic 
oscillator, the distribution is bounded by a circle of radius N/2. 
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